Probing physical properties of confined fluids within individual nanobubbles 
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Spatially resolved electron energy-loss spectroscopy (EELS) in a scanning transmission electron 
microscope (STEM) has been used to investigate as fluidic phase in nanoubbles embedded in a 
metallic PdyoPtio matrix. Using the Is — > 2p excitation of the He atoms, maps of the He distri- 
bution, in particular of its density an pressure in bubbles of different diameter have been realized, 
thus providing an indication of the involved bubble formation mechanism. However, the short-range 
Pauli repulsion mechanism between electrons on neighboring atoms seems insufficient to interpret 
minute variations of the local local measurements performed at the interface between the metal 
and the He bubble. Simulations relying on the continuum dielectric model have show that these 
deviations could be interpreted as an interfzce polarization effect on the He atomic transition, which 
should be accounted for when measuring the densities within the smaller bubbles. 



Confined fluids in nanosized volumes constitute chal- 
lenging objects for both basic and technological aspects. 
The investigation of the structural features and dynam- 
ics of nanojets has given rise to spectacular experimental 
studies and theoretical simulations [lj . Another ideal sys- 
tem is represented by gas confined in nanocavities. It is 
the case of inert gas atoms coalescing as a fluid or a solid 
to fill nanocavities in metals, with spherical or faceted 
morphologies depending of the local pressure. In the 
case of Xe in Al, an interfacial ordering has been demon- 
strated by high resolution electron microscopy [1] . These 
small gas-filled cavities therefore behave as high-pressure 
cells, providing the boundary conditions for the evalua- 
tion of the physical properties of encapsulated gases. A 
most challenging problem is the evaluation of gas density 
and pressure in such cavities. 

Among the possible systems, He nanobubbles in met- 
als have attracted the attention of many researchers, be- 
cause of their high technological interest in the aging of 
the mechanical properties of materials involved in nu- 
clear reactors Q. Measurements averaging the informa- 
tion over large populations of bubbles, the size distribu- 
tion of which being controlled by TEM, have first been 
performed by NMR [4j and by a combination of optical 
absorption and electron energy-loss spectroscopy (EELS) 
0. The first of these studies has revealed a solid-fluid 
transition at 250K for bubble pressures ranging from 6 
to 11 GPa (i.e. He atomic densities from about 100 to 
200 nm" 3 ). The second study comparing UV absorption 
spectroscopy on a synchrotron and high energy resolu- 
tion EELS without spatial resolution on He + implanted 
Al thin foils, have identified the blue shift of the He 
Is — > 2p transition (with respect to its value of 21.218 
eV for the free atom) as a hint for evaluating the local 
pressure . Theoretically, Lucas et al. Q have confirmed 
that this blue shift of the He K-line should be attributed 
to the short-range Pauli repulsion between the electrons 



of neighboring He atoms. Consequently, this effect should 
increase linearly with the density of the He atoms in the 
high-pressure fluid phase likely to exist in these nanosized 
bubbles. Jager et al. 0] have confirmed this linear rela- 
tion between the measured energy shifts (A.E) and the 
average bubble radii (r), the larger shift corresponding 
to the higher He density and consequently to the smaller 
radii. 

With the development of scanning transmission elec- 
tron microscopy (STEM) techniques, capable of measur- 
ing spatially resolved EELS spectra for different positions 
of a sub-nm probe on the specimen, new possibilities 
were offered to perform analysis on individual nanobub- 
bles The most comprehensive study up-to-date has 
been conducted by Walsh et al. Q , who proposed a pro- 
cedure for estimating directly the helium density in a 
single nanobubble. However, this work did not take into 
account the influence of interface excitations on the es- 
timation of the internal density. A more fundamental 
issue which has not been addressed in the case of he- 
lium bubbles is the potential occurrence of density inho- 
mogeneities close do the surface, due to the interaction 
between the confined fluid and the matrix. The investi- 
gation of such effects requires a refined characterization 
at a sub-nanometer scale. 

In this letter, we present a study of the physical param- 
eters (density, pressure, energy of the He K-line) defining 
the state of He inside nanobubbles, by using spatially re- 
solved EELS to map their variations at the nanometer 
scale. The variations between bubbles of different size 
are in agreement with the standard interpretation in the 
literature, while a refined description is required for the 
evolution of the He signal within an individual bubble. 
By using the continuum dielectric model, we show that 
the discrepancies can be explained invoking an effect of 
surface polarization at the interface between the He and 
the metallic surface. This leads to the necessity of a cor- 
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FIG. 1: A: EELS spectra acquired at the centre of bubbles 
of different size (thick lines) and corresponding fit of the Pd 
plasmon (thin lines). B: Subtracted He signal. The shift of 
the He K-line for bubbles of different size is obvious. 

rection of the EELS estimation of the He density inside 
small bubbles. 

The results are issued from an 8-month aged tritiated 
PdgoPtio alloy which exhibits a largely dispersed popu- 
lation of voids (from 2 to 20 nanometers in diameter). 
The EELS measurements have been performed in a VG 
STEM HB 501 with a field emission gun operated at 100 
kV and a homemade detection system formed by a Gatan 
666 PEELS spectrometer optically coupled with a CCD 
camera. Spectrum-images made typically of 64x64 spec- 
tra could then be acquired with the following conditions 
: acquisition time of 200 ms per spectrum, probe of 0.7 
nm with step increments of ranging from 1.5 to 0.5 nm. 
Fig. [TJ^ shows three EELS spectra corresponding to a 
selection of pixels at the centre of three bubbles of dif- 
ferent sizes (Bl, B5 and B7) visible on fig. [2JA These 
spectra correspond to positions where the electron beam 
has crossed both the metallic matrix and the bubbles. 
They all exhibit four major peaks (around 7, 17, 26 and 
33 eV respectively), which are attributed to the low en- 
ergy loss spectrum of the Pd alloy matrix. The sharper 
line between 22 and 23 eV is the signature of the He K 
edge. 

In order to be more quantitative, the He signal of each 
spectrum has been isolated by fitting the palladium con- 
tribution with 4 Gaussian curves (fig. QJ3) . We can then 
identify any change in position, width, total intensity and 
possible occurrence of fine structures or satellites on the 
He K-line, related the different bubbles. For each probe 
position, the He K-line intensity can be evaluated by in- 




FIG. 2: Maps extracted from a spectrum-image of a selected 
area of the sample. A: Bright field image of the analyzed 
area. Bubbles showing He signal are evidenced. B: Helium 
chemical map. C: Map of the He density inside the He filled 
bubbles. D: Map of the energy shift of the He K-line. The 
reference energy is chosen as that of the atomic He (21.218 
eV) 

tegrating the signal over a window of typically 4 eV and 
the results (I He) are displayed as a 2D map (fig. 03) 
of the localization of He atoms. It must be noticed that 
not all the voids contain He atoms. The next step is to 
transform the He K-line intensity map into a cartogra- 
phy of the absolute estimated He density n (expressed 
in atoms/nm 3 ). This can be calculated from the relation 
Q: n = iHe/i&Helzd), where one is the cross section of 
the helium Is — > 2p transition for the used experimen- 
tal conditions (see[9| for calculation), I z is the integrated 
intensity of the elastic peak, d is the local thickness at 
the pixel position of the analyzed He nano-volume. This 
parameter is the source of highest uncertainties. We have 
tested several approaches, but finally we estimated it ex- 
perimentally as the complement to local thickness mea- 
surements of the matrix. The resulting density map is 
shown on fig. 

The mean helium density inside a bubble is estimated 
by averaging the calculated values over a selection of 
pixels corresponding to central positions. The results 
range from 15 to 35 He atoms per nm 3 , the highest value 
been obtained for one of the smaller bubble Bl. The 
energy shift, defined as the difference between the mea- 
sured peak position inside the bubbles and the nominal 
K-line of atomic He @, is mapped on fig. [2p and varies 
from about 1 up to 2 eV. In order to verify the predicted 
linear dependence of AE(n), we have plotted in fig. 
our results issued from several spectrum-images (empty 
squares). A satisfactory fit to a law AE = C n n + D can 
be obtained with C n = (44 ± 7) • 10~ 3 eV • nm 3 and D 
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FIG. 3: A: Experimental relation between energy shift and 
measured density. Empty squares represent uncorrected den- 
sity value; linear fit law: AE = (0.044±0.007)n+(0.07±0.f8). 
Filled square represent density value corrected by surface ef- 
fects; linear fit law:A£ = (0.037 ± 0.004)n + (0.08 ± 0.10). 
Error bars correspond the standard deviation calculated on 
the selection of pixels of the density map, and therefore are 
large for small bubbles having a bad statistics. B: Experimen- 
tal relation between pressure and the inverse of the bubble 
radius. Empty circles are deduced from uncorrected density 
values; linear fit law P = (0.6 ± 0.2) x 1/R ± (0.04 ± 0.05). 
Filled circles account for surface effects; linear fit law P = 
(1.0 ±0.4) x 1/R + (0.002 ±0.083). The theoretical linear re- 
lation for elastic deformation of the PdgoPtio matrix is also 
displayed. Error bars are estimated by calculating the varia- 
tion of the equation of state in the density range defined by 
the corresponding density error bars. Error bars for corrected 
values (not shown) are identical to those for uncorrected val- 
ues 



= 0.07±0.18eK The value of C n lies significantly higher 
than those measured by Jager et al. [3] and Walsh et al. 
[§] but is close to that determined by McGibbon 

Another relevant parameter is the internal pressure. 
In fact, if the bubble deforms the matrix elastically, the 
radius dependence of the bubble pressure is supposed to 
obey an inverse proportionality law P = 2"//r (where 7 
is the surface energy). Following the procedure indicated 
in we calculate the pressure from the measured n by 
using a semi-empirical equation of states (see supplemen- 
tary materials and 

Qmk)- The results are shown in fig. 
[3j3 (empty circles). The pressure inside the bubbles is 
shown to increase roughly from 0.1 to 0.3 GPa (i.e. in a 
range well below the solid to liquid transition pressure), 
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FIG. 4: Experimental profiles of the estimated density ( filled 
circles) and of the blue-shift of the He K-line (filled squares) 
as a function of the mean local thickness d. Data are extracted 
from a spectrum image of 40x40 pixels (spatial sampling of 
0.5 nm), acquired on a bubble of 19.5 nm diameter For com- 
parison, the corresponding simulated profiles of the He-K line 
density (empty circles) and energy position (empty square) are 
also displayed. 



when the diameter of the bubble decreases from 17 to 
5 nm. A reasonable value for the surface energy of the 
PdcjoPtio alloy is 7 = 1.9Jm -2 , to be compared to our 
experimental slope 0.3Jm~ 2 . Then, the bubbles seem to 
be under-pressured at the moment of our TEM observa- 
tion. 

The spectrum-image technique offers the extra possi- 
bility of exploring any potential intra-inclusion spatial 
dependence. A varying contrast seems visible within 
larger bubbles in fig. [2jZ3 and [2jD, where the density 
drops while the energy shift increases close to the sur- 
face of the bubbles. In order to further investigate this 
behavior, fig. H] shows experimental profiles of density 
and shift as elaborated from a spectrum-image of a 19.5 
nm bubble, probed with a better lateral sampling of 0.5 
nm. Each point of the two profiles has been calculated 
selecting annular regions of pixels corresponding to the 
same analyzed thickness, and fitting the He Is — > 2p 
transition with a Gauss function to calculate the energy 
position and the intensity. From the centre to the bub- 
ble surface a 37 percent drop is observed for n while the 
energy shift increases by 0.17 eV, which is one order of 
magnitude smaller than the shift between different bub- 
bles. This anticorrelation is in contradiction with the 
general tendency previously observed between individual 
bubbles. Indeed, when only Pauli repulsion between He 
atoms is taken into account, such a density drop should 
lead to a 0.35 eV shift toward lower energies. In order 
to evaluate the potential occurrence of surface effects, we 
have performed EELS spectra simulations, by adapting 
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to the case of embedded spheres the continuum dielectric 
model which has proved its efficiency for modeling local 
surface phenomena in nanosized systems, such as single- 
walled nanotubes 11, ill- As an input for the simulation, 



we used a lorentzian dielectric constant corresponding to 
a He fluid of constant density Q. After simulation of 
spectra for different local thicknesses d, the procedure 
used to extract n and AE on experimental data is ap- 
plied. The resulting simulated profiles are compared to 
the experimental ones in fig. [4)3 . Both AE and n varia- 
tions are reproduced but underestimated. Consequently, 
the major part of the effect can not be attributed to 
a real change in the density, since the model assumes 
a constant one, but to the influence of surface excita- 
tions on the measurement. We stress that the evidenced 
surface effect is not due to an usual plasmon mode, be- 
cause it does not correspond to a pole (resonance) but 
to a maximum of the dielectric response of the sphere 
Im{{e He -e m ) / {e He +2e m )) (where e He and e m are the di- 
electric constants of He and of the metallic matrix respec- 
tively). An interface plasmon excitation is expected at a 
lower energy value (of the order of 7 eV), and is rather 
unsensitive to the helium density. However, the dielectric 
formalism commonly used to model plasmon excitations 
furnishes reliable (similar) interpretations for the effects 
of interface polarization on the atomic transition. It is 
well known that such a formalism is very sensitive to the 
input dielectric constants of both materials, and discrep- 
ancies between simulated and experimental data can be 
partially explained by a lack accuracy of the lorentzian 
model adopted for He as well as of the experimental data 
used for Pd. Nevertheless, the energy shift of He K-line is 
explained by the contribution of a surface "mode" which 
energy, for this particular system, is slightly higher than 
that of the bulk He line (see supplementary materials). 
The decrease of the estimated density can be related to 
a companion effect known, for plasmons, as a boundary 
effect or "Begrenzung" effect [13j. This effect is com- 
monly interpreted as a modification in the probability to 
excite bulk modes due to the occurence of surface excita- 
tion, and reveals itself as a negative contribution to the 
intensity of the bulk He line, the importance of which in- 
creases as the He thickness decreases. We point out that 
this is the first time that this surface effect, which has 
been thoroughly investigated for valence electron excita- 
tions, is evidenced on an atomic-type excitation, using 
EELS. 

Therefore, beside its intrinsic fundamental interest, 
this surface-induced decrease in the density estimated 
from atomic transition signal should be taken into ac- 
count in the study of the bubble formation mechanism. 
We calculated a correction coefficient G to apply to ex- 
perimental intensities in order to account for surface ef- 
fects in the estimation of n. Such a coefficient is given 
by the ratio G = I WO s/Itot, where I wos is the He-K in- 
tensity simulated excluding surface contributions (ideal 



case), and Itot is the total simulated intensity (real case). 
The resulting AE(n) corrected relation is displayed in 
fig. [3K (filled squares). The linear fit gives an estima- 
tion of the slope C„ decreased of 19 percent and closer 
to the values in the literature. Even larger is the cor- 
rection to relation between the pressure and the inverse 
radius (filled circles in fig. [3)3), with a slope increased 
by a factor close to 2. Nevertheless the comparison of 
the corrected data-set to the linear relation characteris- 
tic of the elastic deformation regime (also displayed in 
[5)3) confirms that the bubbles are under pressured. 

In conclusion, the present analysis of the confined He 
fluidic phase, at their interface with the embedding ma- 
terial, has evidenced an interface-induced effect on the 
atomic-like spectral transition in He. Consequently, a re- 
liable estimation of the helium internal density and pres- 
sure requires a correction from this surface effect, espe- 
cially in the case of small bubbles. This effect can be 
of much broader interest, and it should also be identi- 
fied in quite different situations, such as those encoun- 
tered on semi-core loss edges (HLO2.3) in dielectric thin 
films [3] .From another point of view, the interpretation 
of residual discrepancies between experiments and simu- 
lations, when exploring the influence of the distance from 
the interface, require further modeling, accounting for 
changes in Pauli repulsion and in Van der Waals forces 
close to the interfaces. 

We thanks L. Henrard, A. A. Lucas, Ph. Lambin and 
P. Loubeyre for fruitful discussions 
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